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UTILIZING TYPE lA SUPERNOVAE IN A LARGE, FAST, IMAGING SURVEY TO CONSTRAIN DARK ENERGY 
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ABSTRACT 

We study the utility of a large sample of type la supernovae that might be observed in an imaging survey 
that rapidly scans a large fraction of the sky for constraining dark energy. We consider both the information 
contained in the traditional luminosity distance test as well as the spread in la supernova fluxes at fixed redshift 
induced by gravitational lensing. As would be required from an imaging survey, we include a treatment of 
photometric redshift uncertainties in our analysis. Our primary result is that the information contained in 
the mean distance moduli of supernovae la and the dispersion of supernova la distance moduli complement 
each other, breaking a degeneracy between the present dark energy equation of state and its time variation 
without the need for a high-redshift (z}t 0.8) supernova sample. Including lensing information also allows for 
some internal calibration of photometric redshifts. To address photometric redshift uncertainties, we present 
dark energy constraints as a function of the size of an external set of spectroscopically-observed supernovae 
that may be used for redshift calibration, A^spec- Depending upon the details of potentially-available, external 
supernova data sets, we find that an imaging survey can constrain the dark energy equation of state at the epoch 
where it is best constrained Wp, with a Icr error of a(wp) w 0.03-0.09. In addition, the marginal improvement 
in the error a(wp) from an increase in the spectroscopic calibration sample drops once A^spec ~ a few x 10"^. 
This result is important because it is of the order of the size of calibration samples likely to be compiled in the 
coming decade and because, for samples of this size, the spectroscopic and imaging surveys individually place 
comparable constraints on the dark energy equation of state. In all cases, it is best to calibrate photometric 
redshifts with a set of spectroscopically-observed supernovae with relatively more objects at high redshift 
(z^ 0.5) than the parent sample of imaging supernovae. 
Subject headings: 



1. INTRODUCTION 

Type la Supernovae (SNela) can be used as well-calibrated 
standard candles with rel atively little dispersion in their in- 
trinsic luminosities (e.g.. lPhill ips"1993': 'Ilamuv et aLlll995l; 
l^ess et al. 1996; Pri eto et al.l2 006; Guy et al. 20071: iJha et all 
|2007). This has allowed for recent measurements of the rela- 
tion between cosmological distance and redshift, which pro- 
vide the strongest contempora ry evidence for an accelerating 
cosmologic al expansion (see lAstier et al.l 120061: iRiess et akl 
12007; Woo d-Vasev et al.ll2007l for recent applications). Sev- 
eral studies have examined the spread in observed fluxes 
(rather than their average) due to the magnification of SNela 



die distance-redshift test ("Sasaki' 119871; 


Linder et al.l 


1988 


Kantowski et al. 1995; Wambsganss et al. 


19971 Holz) 


1998 



Wangj il999t iBarben i2000i: |HoIz & Linder 2005), a cross- 
check of lensing shear maps (jCporay et al. 2006j and refer- 
ences therein), or as an interesting signal in its own right 
jMetcalf '1999^, 'Dodels on & VaUinottol l2006t ICoorav et all 
12006; Albrecht et al. 200$ 

The present paper is inspired by these studies and other ef- 
forts to explore how SNela might be utilize d as part of a large, 
photometric survey (e .g., Albrec ht et al.l 12 006; Zha n et al.l 
120081; iHannestad et"ai] 12008; .Zhang & Chen .2008.) . The 



aforementioned studies considered the lensing of a relatively 
small (;$ 10^) spectroscopic sample of SNela, with redshifts 

tha t are known p recisely (with errors ^ 10"^, as prescribed 
bv lHuterer et alJIIOOl . We examine the utility of a signifi- 
cantly larger (~ 10^) photometric sample of supernovae that 
may be collected by a future, wide, fast, imaging survey such 



as the proposed Large Synoptic Survey Telescope (LSST)' in 
which the vast majority of SNela will not have accompanying 
spectroscopic observations. It is thought that the identification 
and determination of ph otometric redshif t s for SNela can be 
done relatively reliably jPinto et al.ll2004i; IWang et al.| [2007). 
Consequently, the dispersion in SNela fluxes at fixed photo- 
metric redshift has three sources: the intrinsic dispersion of 
SNela, including errors in calibrating the standard candles 
and applying K-corrections and extinction corrections; pho- 
tometric redshift errors; and any dispersion induced by grav- 
itational lensing. Moreover, each of these factors are thought 
to contribute to the total dispersion at fixed photometric red- 
shift at comparable levels (see §|2l). The implication is that a 
very large sample of supernovae from an imaging survey can 
be used to overwhelm the intrinsic dispersion, thus enabling 
the measurement of SNela dispersion induced by photomet- 
ric redshift errors and gravitational lensing which should shed 
light on both SNela photometric redshifts and cosmological 
parameters at relatively low redshift. 

The modest goal of this paper is to make preliminary es- 
timates of the utility of such a photometric SNela sample 
to learn about cosmology from the dispersion among SNela 
fluxes in conjunction with the traditional distance-redshift 
test. We find that the spread in observed SNela fluxes at 
fixed redshift provides constraints that are not, by themselves, 
competitive with other probes (such as weak lensing shear), 
but are yet interesting, subject to unique systematic errors, 
and obtained without significant additional observational ef- 
fort if such a sample is to be used to perform the traditional 
distance-redshift test or to map baryon acoustic oscillations 

' |http : / /www. Isst ■ orgi 
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(as studied in lAlbrecht etani2006l: Izhan et al.ll2008l) . In fact, 
because such a measurement will be most constraining on 
dark energy at low redshift (z ~ 0.2-0.4) it will have some de- 
gree of complementarity with other measures, such as cosmic 
shear which best constrains the dark energy equation of state 
near z ^ 0.7. More importantly, the information contained 
in the traditional distance-redshift test (most sensitive to the 
dark energy equation of state near z ~ 0.1) and the disper- 
sion of SNela brightnesses complement each other and help 
to make cosmology with photometric SNela more robust to 
uncertainties in photometric redshifts. 

Of course, the lack of spectroscopy is a major complication 
in using SNela from a photometric survey. As the use of pho- 
tometric SNela samples for cosmology requires some treat- 
ment of photometric redshift uncertainties, we also address 
the dependence of cosmological constraints upon the calibra- 
tion of SNela photometric redshifts. In particular, we find 
that utilizing the dispersion in SNela distance moduli allows 
for mild self-calibration in the uncertainties of photometric 
redshift distributions. Perhaps more interestingly, we present 
results for dark energy constraints as a function of the size of 
a spectroscopic sample of SNela that can be used to calibrate 
photometric redshifts for two different assumptions about the 
redshift distribution of this spectroscopic calibration set. This 
is a convenient way to express cosmological constraints as a 
function of prior knowledge of the photometric redshift dis- 
tribution of SNela. Interestingly, we find that if the goal is 
to seek the best constraint on the equation of state parameter 
at the epoch where it is best constrained by the data (the so- 
called "pivot" equation of state, see §|2] results are similar for 
a constant equation of state model), then there is decreasing 
marginal value in expanding the size of the spectroscopic cal- 
ibration set beyond a few thousand SNela, a number of SNela 
that should be achievable in the coming decade. However, 
this statement does depend upon model details and upon the 
assumption that it is the pivot equation of state that is of most 
interest. 

There are several important caveat to our analysis. We 
neglect the, perhaps considerable, additional complications 
in supernova type identification and calibration arising from 
the lack of spectroscopy. Aside from neglect of several 
potentially-important observational realities, our calculation 
is also an idealized one. We have treated the influence of grav- 
itational lensing on observed fluxes in a simplified manner, in 
large part because current analytic models for the proposed 
lensing signal are not adequate for future surveys and a full 
treatment of this effect is computationally demanding. How- 
ever, part of the purpose of this paper is to demonstrate that 
such a signal does contain valuable information and to empha- 
size that as we yet lack the theoretical tools needed to perform 
a robust analysis of future data, developing such tools should 
be a high priority in the run-up to forthcoming photometric 
surveys. 

The present paper is organized as follows. In § |2]we de- 
scribe the our treatment of SNela observables and photomet- 
ric redshifts as well as our methods for forecasting cosmolog- 
ical parameter constraints. We present our primary results in 
§|3] We discuss our results, including implications and impor- 
tant caveats in § |4] 

2. BASIS AND METHODS 

The average distance moduli /i, of high-redshift supernovae 
as a function of redshift provide the most direct evidence for 
cosmological acceleration. In the absence of any lensing, the 



distance modulus of a standard candle at true redshift z is 
^ = 5 log[DL(z)/10 pc], where Z)l(z) is the luminosity distance 
to redshift z- We explore both the mean and dispersion in dis- 
tance moduli of SNela measured in a future, large photometric 
survey. We assume a survey that covers w 20,0000 deg- to a 
redshift of z ^ 0.8 for 10 years as might be achieved with the 
LSST instrument and surve y strategy. We tak e a SNela red- 
shift distribution dn / dz as in lzhanetal.1 (l2008h . which rises to 
z ^ 0.5 (dn/dz cx z^ ' for z <C 0.5) and declines rapidly there- 
after (dn/dz oc exp[-32(z- 0.5)2] f^j. ^ ^ q 5-, n^^^ ^^^^ 
which supernovae may be discovered and the efficiency with 
which they will produce light curves suitable for calibration 
and reliable photometric redshift determination is uncertain 
and depends upon the exposure strategy of any such survey. 
For simplicity, we choose a fiducial survey containing one 
million SNela. We note in passing that this normalization 
is significantly more ambitious than the 3 x lO'' SNela as- 
sumed for the hypothetical, stage IV, ground-based survey ex- 
plored by the Dark Energy Task Force (DETF, Albrecht et al] 
|2006j), but significantly less ambitious than the 7.4 million 
SNela assumed in the explor ation of SNela as a tracer of 
baryon acoustic oscillations bv lZhan et al.l (l2008h . In the ab- 
sence of additional prior constraints on cosmological param- 
eters and external calibration samples of SNela (see details 
below) parameter uncertainties would scale with sample size 
as CT oc 1 / V-^SNe; however, our constraints scale more slowly 
because we assume priors on cosmological parameters as well 
as the availability of additional spectroscopic SNela samples 
that improve the calibration of SNela properties and redshifts. 

Each observed supernova will be assigned a photometric 
redshift Zp that is significantly less reliable than could be 
achieved with spectrosco py. We model photometric redshift 
uncertainties as follows ( Ma et al.1l2006l) . We assume a distri- 
bution of photometric redshifts given a spectroscopic redshift 
P{z^\z), that is a Gaussian with an offset from the true red- 
shift b^{z) and dispersion (t,-(z), each of which may depend 
upon redshift. Recent studies indicate that b- ^ and a? ~ 
(T.o(l +z) with (7,0 Ki 0.01 sho uld be achieva ble (iPinto et alJ 
l2004t lAlbrecht et all bOOe; Wa ng et al.ll20O7l) . In practice, it 
is necessary to allow and ct; to have significant freedom to 
vary as a function of z. We do this by representing the func- 
tions br and cr, by the values of these functions tabulated at «pz 
redshifts spaced evenly in z from z = to z = 1 .5 and evalu- 
ating these functions at points between the redshifts at which 
they are tabulated using linear interpolation. Thus there are 
2npz photometric redshift parameters. We set «pz by requir- 
ing that Hpz be large enough that when all 2npz parameters 
are allowed to vary without prior knowledge, the constraining 
power of the data converges to a minimum value correspond- 
ing to complete lack of redshift information. In practice, our 
results converge to this minimum with «pz w 15 and t o be con- 
servati ve, we fix Mpz = 3 1 as in the lensing study of iMa et alJ 
( I2OO6I) for the remainder of this paper 

Of course, there will be some prior knowledge of f (Zp |z) 
and in principle each of the photometric redshift model pa- 
rameters will have its own prior In fact, some calibration of 
the photometric redshift distribution is needed in order to uti- 
lize SNela from an imaging survey as we demonstrate explic- 
itly in the following section. For simplicity, we quantify the 
influence of priors using two one-parameter models that are 
reasonable, but not exhaustive of all possibilities. We assume 
a sample of A^spec SNela with known spectroscopic redshifts. 
We further assume that this calibration sample either follows 
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Fig. 1 . — The dispersion in supernova distance moduli. Left: Dispersion as a function of photometric redshift for SNela in 51 bins of photometric redshift. The 
horizontal, dotted hne represents the assumed intrinsic dispersion level of (T|„( =0.1. A set of three, nearly indistinguishable dashed lines decrease as a function 
of redshift. These lines represent the contribution to the dispersion from the first two terms in Eq. (4). The central line corresponds to our fiducial cosmological 
model while the upper (lower) line corresponds to the same cosmology but with vvq = — 1.1 (vi'o = —0.9). A set of three solid lines increases monotonically with 
redshift. These represent the contribution to the net dispersion from gravitational lensing. The central line is for our fiducial cosmology, while the upper (lower) 
line has wq = — 1 ■ 1 (wo = —0.9). The uppermost, solid lines represent the total dispersion in Eq. (4)- The central, upper, and lower fines correspond to models with 
Wo = -1, Wo = -1.1, and wq = —0.9 respectively, while the errorbars correspond to the statistical error on a determination of each of the 51 cf^\ from a sample 
of one million SNela. Right: Residuals of the dispersion in distance modulus with respect to our fiducial cosmological model with wg = — 1. The upper dashed 
line is the dispersion in a cosmology with »'o = -1.1 and the lower line is the dispersion with wq = -0.9. The e rrorbars are the en'ors on the measurement of the 
dispersion in each of 51 photometric redshift bins as in the left panel. This figure is modeled after Figure 1 of lDodelson & Vallinotto 12006) and demonstrates 
that the dispersion in standard candle distance moduli that may be measured in a forthcoming, large, photometric survey like the LSST will undertake, has some 
leverage to distinguish dark energy equations of state. 



the overall distribution An /dz, or is spaced uniformly in red- 
shift from z = 0. 1 to z = 1 .7 as might be expected from a Joint 
Dark Energy Mission (IDEM) SNela sample^ that might be a 
conte r nporary of a large-sca le photometric survey (iKim et alj 
l2004t lAlbrecht et"aLl (2006h . We also incorporate the influ- 
ence of a nearby sample of SNela with spectroscopic redshifts 
that help to serve as a low-redshift anchor for the distance- 
redshift test. For this, we take the nearby SNela sample to 
contain 500 objects uniformly distributed in the redshift inter- 
val z = 0.03-0.08 with redshift errors of a- = 1 0"'* as expected 
from the ongoing Nearby Supernova Factory (lAldering et al.l 
|2b02; Copin et al. 2006; Albrecht et al. 2006)^ experiment. 

Any calibration set would determine the bias, b\, and dis- 
persion, of the photometric redshift distribution in red- 
shift bin i with uncertainties of cr^i = a\/ y^A^j^ and (t„, = 

cr^/y^2A^^pg^ where A^pec is the number of spectroscopic SNela 

that fall in bin i. We include the influence of prior knowledge 
of the photometric redshift distribution using such priors. We 
assume a fiducial model set by b, = and a, = <J:o(l +z) with 
iT;0 = 0.01. We reiterate that, unless otherwise stated, we al- 
low for uncertainty in the photometric redshift distribution as 
a function of redshift through the 2npz photometric redshift 
parameters and do not assume that the photometric redshift 
paramete rs are fixed at the fiduc ial values (in contrast to, for 
example, iHannestad et al.ll2008l who assum e perfect knowl- 
edge of SNela redshifts and lZhan et al.ll2068l who do not con- 

^ The Supernova Acceleration Probe (SNAP) is an example, see 
|http : / / snap . Ibl . g ov| 

-' Ihttp : //snf actory . Ibl ■ govl 



sider calibration of the bias in photometric redshifts). 

For computational convenience, we compute the observ- 
able properties of supernovae in bins of photometric redshift. 
Given the distribution P(zp|z), the true redshift distribution of 
SNela in photometric redshift bin i is 

dn\ f'f dn 

where zjf™ and Zp'°'' are the lower and upper boundaries of 
photometric redshift bin i respectively. The total number of 
SNela in photometric redshift bin i is N'^^^^^ = /(dni/dz)dz and 
we define the normalized redshift distribution of SNela in bin 
i as gi(z) = (d«i/dz)/A^phof In practice, we take 51 bins spaced 
equally in redshift from z = to z = 1.5. Though convenient 
considering observables averaged within redshift bins reduces 
information content relative to the total information available. 
We have explored various binning schemes and find that finer 
binning is unnecessary as our results are already converged 
with this number of bins. In particular, increasing the number 
of bins by a factor of 3 gives only ^ 4% better constraints 
on the pivot dark energy equation of state parameter defined 
below. 

Consider an effective distance modulus defined for conve- 
nience to be the difference between apparent magnitude and 
an assumed average absolute magnitude. Any observed su- 
pernova has an effective distance modulus 

fi = A + 5log( -/^]+Sfii„i + dfiisn>i, (2) 
VlOpcy 

where (5/iint is the variation in apparent magnitude due to the 
intrinsic variations in SNela absolute magnitudes, Sfij^ns is the 



4 



Zentner & Bhattacharya 



contribution from gravitational lensing, and A is some con- 
stant offset that is applied to all supernovae and represents 
a deviation from the assumed average absolute magnitude 
which includes errors in Hq. In what follows, we set A to 
vary freely. The average distance modulus of SNela in the i* 
photometric redshift bin is 



Ati = A + 5 / dzgiiz)\og 



10 pc 



(3) 



and the standard distance-redshift test is an application of this 
relation. In addition to distance modulus, the Greek letter "/i" 
is also the conventional symbol for the lensing magnification. 
In order to avoid confusion, we will not refer to the lensing 
magnification directly but only to the shift in distance modu- 
lus induced by weak gravitational lensing. 

The dispersion in the effective distance modulus cr^ i, is 
given by 



<i = 25 



dzgiiz)\og 



dzgiiz)\og 



D, 



10 pc 



10 pc 

+ / dzgi(z)al,(z)+ I dz giiz)crl„,iz) 



where aint is the intrinsic dispersion of SNela brightnesses, 
which may itself vary with redshift, and ciens is the disper- 
sion induced by gravitational lensing. The first two terms in 
parenthesis in Eq. (01 represent the spread in observed SNela 
distance moduli due to the spread in photometric redshifts. 

In the weak lensing limit (c onvergence k <^ 1), the disper- 
sion due to lensing is ( Bernard eau et alJI 19971 : IValageasll2000l : 
iDodelson & VallinottQ,2006.) 



2 

"lens 



iz) = 



87rW(10) 



dz 



, WHz',z) 
H(z) 



dkkP(k,z')W^HiDf,ke,), 



(5) 



where W(z',z) = HoDa(z')Da(z\z)/Da(z), Da(z) is the angu- 
lar diameter distance to redshift z, and Da{z',z) is the angu- 
lar diameter distance between redshifts z' and z- The quan- 
tity Pik,z) is the matte r power spe c trum and we evaluate 
it using the relation of ISmith et all (l2003l) . The function 
yVjuix) = 2Ji(x)/x in Eq. (|5]l is a smoothing function that 
arises by considering the magnification averaged over an an- 
gular tophat of radius 1. In practice, we set = 1", but 
the choice is relatively unimportant in that DAk9^ ^ 1 and so 
y^miDAkO^) ~ 1 for relev ant choices of 9^. We assume the 
lensing to be uncorrelated (lMetcailll999l:ICooray et al.ll2006l: 
IDodelson & Vallinottoll2006l) ' and take the tr^ i as an indepen- 
dent set of observables that can be extracted from the observed 
SNela population. 

Figure [T] shows the dispersion in distance modulus in 
51 photometric redshift bins as a function of the center of 
the photometric redshift bin z-p in our fiducial cosmologi- 
cal model. Our fiducial model has t^M = 0.13, wb = 0.0223, 
Vlk = 0, A^ = 2.47 X 10"'', and n., = 1 . We describe dark energy 
by its present energy density i^DE = 0.76 and a time-varying 
equation of state parameter w{a) = - a), with wq = -\ 

and Wa = in our fiducial model. For the purpose of this illus- 
tration, we take (Tint = 0.1. We will return to the intrinsic dis- 
persion later Notice that the contributions to the dispersion 
due to the photometric redshift distribution and gravitational 
lensing are comparable and give the total dispersion a charac- 
teristic dependence upon redshift. In particular, the dispersion 



decreases at low redshift and increases at high redshift due to 
the increased dispersion due to lensing. This makes it possi- 
ble to use the dispersion to extract cosmological information 
without a good constraint on the level of the intrinsic disper- 
sion. 

We estimate constraining power using a Fisher matrix 
analysis. The Fisher matrix formalism is ubiquitous in 
cosmological parameter forecasting (useful referenc es in- 
clude [Jungman et al. 1996; Tegmarket al. 1997; Seliakl l997t 
Metcalf 1999t lAlbrecht et al.ll2006tlZhan et al.ll200 8. the last 
three of which explore cases similar to the present paper) and 
so we will not review the formaUsm here. We take as our ob- 



servables the set /ii and cr^ j with errors of a{p.i) = cr^.i/ ^N^piioi 

and (T((T^i) = cr^j/ ^JlN^^ respectively. We refer to con- 
straints using the p.^ alone as "luminosity distance test," and 
constraints from the cr^ i alone as the "dispersion test." The 
2^ Fisher matrix F,„„, approximates the covariance in model pa- 
rameters locally about the maximum likelihood, so the in- 
dices run over the parameters in the parameter space we seek 
to constrain. The la constraint on parameter p is approxi- 
(4) mated by the square root of the diagonal component of the 
inverse Fisher matrix corresponding to this parameter We 
denote this <T{p) ~ y/[F~^]pp. Along with the parameters 
A and aint describing SNela and the parameters describing 
the photometric redshift distribution, we vary the eight cos- 
mological parameters in the previous paragraph about their 
quoted fiducial values. We take independent, Gaussian pri- 
ors on these parameters of a{ujM) = 0.004, (t(cjb) = 6 x 10""*, 
airtk) = 10"^ cr(ln A^) = 0.04, and a(n,) = 0.02, all of which 
are c omparable to contemporary constraints on these parame- 
ters (iKomatsu et al. 2008). 

Absent significant direction from theory, it is not clear 
how to assess the constraining power of a dark energy pro- 
gram. We use two reasonable metrics. First, we consider 
the area of the marginalized 95% error ellipse in the wo-vva 
plane A, as suggested by the DETF ( Albrecht et al. 2006L 
the DETF actually uses 1 /A, see also Huterer & Turner 200T| 
for a similar suggestion). Alternatively, it may be supposed 
that the goal of a dark energy experiment should be pri- 
marily to limit any deviation from a vacuum energy or cos- 
mological constant, both of which have constant w = -I. 
In this case, it is interesting to examine the error on the 
equation of state parameter w{a) at the epoch where it is 
most well constrained, the so-calle d pivot scale factor a„ or 
pivot re dshift Zni y = 1 /flp - 1 (f g-. iHuterer & Turned 120011 : 
iHu & J ain 2004; lAlbrecht et al]|2006h . The pivot scale fac- 
tor is Qp = 1 + [i^"']w„w,/[f'~^]waH'a- Thc corrcspondlng pivot 
equation of state is Wp = wo + (l-ap)Wa, so the error on Wp is 

cr^(wp) = [F'^lvowo - [^"']^„»va/[f'"^]w,wfl- The parameters Wp 
and Wa are uncorrelated and the transformation from a wq-w^ 
parameterization preserves the area of the error ellipse in the 
two parameters, so ^ ~ 6.177r(T(Wp)cr(wa). 

Our estimate of parameter uncertainties stands as one of 
many abuses of the Fisher matrix that appear in the litera- 
ture and this warrants some discussion. For one thing, we 
expect the likelihood to be non-G aussian for an interesting 
range of parameter values (though Metcalf 1999', also works 
in the Gaussian approximation, see also § |4|. More impor- 
tantly, we consider only the mean [Eq. ©I and dispersion 
[Eq. ©I, rather than the full /i distribution. In principle, one 
would like to, and should, account for the full shape of the 
distribution of /i, which is non-Gaussian in large part because 




Fig. 2. — Constraint contours in the wq-w^ plane. The left panel shows la constraint contours in the limit of perfect knowledge of the parameters of the 
photometric redshift model. The right panel shows parameter constraints assuming priors on the photometric redshift parameters that would be obtained after 
calibrating to a spectroscopic sample of A'spcc = 3000 SNela distributed uniformly in redshift over the interval 0.1 < z < 1.7. Note that the axes in the two panels 
span different ranges. In each panel, the cross marks the fiducial model with wq = —1 and = 0. The outer, doshed contours represent the dark energy constraints 
from the standard luminosity distance-redshift test [Eq. (3j] alone. The outer, dotted contours represent constraints from the dispersion among SNela alone. The 
solid contours show the combined constraints from both the mean and dispersion among SNela distance moduH. To set these constraints in context, the inner, 
thin, dotted contours in each panel show the constr aints expected from galaxy weak lensing shear tomography in an LSST-like survey computed using the same 
fiducial model as described in lZentner et al.l 420081) . 



the distribution of magnifications due to lensing is strongly 
non-Gaussian (e.g., Schneider & Wagoner 1987; Sasaki 1987; 
Wambsganss et al. 1997; Holz 1998; Wang 1 19991; IValageasI 
2000; .Munshi & Jain 2000; .Wang et al.,.2002h . One way to 
see that lensing should be non-Gaussian is to recognize that a 
minimum amount of de-magnification (or dimming) will oc- 
cur when the null geodesic passes through an empty beam 
with density p = 0. The minimum magnification of a source at 
redshift Zs, corresponds to a maximum shift in distance mod- 
ulus of (5^max = (15/2)riMf^o W {z,Zs) / H{z). There is no 
corresponding upper limit to magnification or corresponding 
lower limit to (5/iiens- 

We use the Fisher matrix formalism and consider only 
the dispersion in /i as a matter of pragmatism. First, com- 
puting the distribution of magnifications due to lensing is 
computationally intensive and still subject to uncertainties in 
numerical mo deling of nonlinear structure formatio n (e.g. , 
iHuterer & Takad a 2005; Hagan et al.ll2005i: iRudd et al.ll20C)8l) 
at levels that are important for forthcoming data. Neglect- 
ing systematic issues, this problem could be circumvented 
if a reliable, analytic fitting form for the magnification dis- 
tribution could be used. Wang et al. (2002) provide such a 
fit, but we find that the IWang et alJ (120021) relation is neither 
sufficiently accurate to address forthcoming large data sets 
(~ percent-level predictions are necessary), nor is it internally 
self-consistent below redshifts z ~ 0.6, where the probabil- 
ity density of magnification exhibits discontinuities and vio- 
lates flux conservation'*. This is particularly important for the 

^ The redshifts at which these issues of inconsistency become important 
vary considerably with cosmological parameters. These issues are alluded 
to in the work of Dodelson & ValUnotto ( 2006), but they do not specify the 
shortcomings of the Wang et al. (2002) fitting form. The shortcomings are 
understandable in the sense that the present application was likely not fore- 



present study because the bulk of SNela in any ground-based, 
photometric survey of a large fraction of the sky will be at 
redshifts < 0.6. 

If such a fitting form were available, the present study 
with 10^ SNela and ~ 70 parameters would still be 
computationally demanding. What one would like to do 
is to perform a Monte-Ca rlo analysis (as was done in 
iDodelson & Vallinottol |20()6|) in which random realizations 
of SNela fluxes are generated and to map out the likeli- 
hood during a random walk through the parameter space 
(e.g., [C hristensen & Mever '2001'; Chr istensen et all 120011; 
iKnox et al. 2001; Kosowskv et al. 2002). In doing so, one 
would also treat the full lensing distribution, rather than tak- 
ing the weak lensing limit (where convergence k ^ 1 or 
(^Miens ^ 1), which does not account properly for the ob- 
jects in the high-magnification tail of the lensing distribution. 
Dodelson & Vahinotto (2006) and Sarkar et al. (2008b) have 
shown that working in this limit does induce biases in cosmo- 
logical parameters inferred from any observational data sets, 
but does not significantly alter parameter constraints. Indeed, 
we note here that we have successfully repeated the Monte 
Carlo analysis of Dodelson & Vallinotto (2006), finding sim- 
ilar results. We have also performed the same Monte Carlo 
analysis for our photometric survey in cases of a reduced pa- 
rameter space where the photometric redshift distribution is 
assumed to be known perfectly and in which only Om and 
W() vary. In both cases, a Fisher matrix analysis yields con- 
straints that are in rough agreement (^ 30%) with the more 
complete approach. Detailed and precise theoretical predic- 
tions and detailed survey strategies and systematics analyses 

seen by IWang et alj j2002f) . so it seems reasonable that these authors would 
not expend significant effort to calibrate their fit at low redshifts. 
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are not yet available, so we do not consider a calculation that 
incorporates a full Monte Carlo analysis of the lensing distri- 
bution and the parameter space to be warranted at present and, 
for practical reasons, we have not performed such a calcula- 
tion. However, any analysis of an observational data set must 
undertake such a calculation. 

Rather than treating such a specific calculation, our point is 
to indicate that the dispersion among supernovae in any large 
photometric survey that also aims at using the SNela sample 
to perform the distance-redshift test, measure baryon acoustic 
oscillations, or even assess SNela systematics will be a source 
of meaningful information. In fact, part of our point is that 
such information may be present in future data, but the theo- 
retical tools to do a rigorous analysis of such data are lacking 
and must be developed. We return to many of these issues in 

m 

3. RESULTS 

This section contains our results regarding the utility of a 
large, photometric sample of SNela to constrain cosmological 
parameters. Our primary interest is in the ability of SNela to 
constrain the dark energy equation of state, and so we focus 
our attention on wo, Wa, the pivot value of the equation of state 
parameter Wp, and the 95% wq-w^ ellipse area A. We illustrate 
the general utility of SNela from a large photometric survey, 
and in particular the utility of examining both the mean and 
dispersion among distance moduli, in Figure |2] 

The left panel of Fig. |2] shows constraints in the 
unrealistically-optimistic situation where the parameters of 
the photometric redshift model are known perfectly. In real- 
ity, this cannot be the case, but this represents the limit of the 
best possible cosmological constraints that could be achieved 
with a photometric survey of supernovae. We present results 
for a fiducial case with lO*" supernovae. We reiterate that 
for a fixed redshift distribution and in the absence of exter- 
nal priors, constraints would scale with the total number of 
supernovae as u oc 1 /V-^SNe and this could be used to scale 
constraints to approximate those from smaller or larger sam- 
ples with similar redshift di s tributions (e.g.. [ Albrecht et alj 
I2006t iHannestad etal] 120081; IZhan et al.l l2008l) . Our c"or> 
straints scale somewhat more slowly due to the influence of 
the external priors we assume. 

There are several things to note in the left panel of Fig. |2l 
First, in the limit of perfect knowledge of the photometric 
redshift distribution and perfect control of systematic uncer- 
tainties, constraints from the standard luminosity distance test 
dominate, placing a constraint on Wp of a{wp) ~ 0.01. How- 
ever, the extra information contained in the SNela dispersions 
is not uninteresting. The dispersion test constrains the pivot 
equation of state parameter to <j(wp) ~ 0.07 by itself. More 
importantly though, is that the luminosity distance and dis- 
persion contain complementary information. The luminosity 
distance test alone has a pivot redshift of z 0.1, and so the 
constraint contours from this test alone are quite vertical. The 
dispersion test has a pivot redshift at z ~ 0.4, and so the disper- 
sion contours are more inclined and it is clear that the disper- 
sion test complements the distance test. The combined con- 
straint from the mean and dispersions among supernovae dis- 
tance moduli result in a constraint on Wp that is only slightly 
smaller than that from the distance test alone, but a significant 
reduction in the area of the 95% contour by a factor of roughly 
~ 2. For reference, we also show constraints from weak lens- 
ing shear tomograph y for the LSST survey as estimated by 
IZentner et al.l (l2008l) . assuming photometric redshift calibra- 



tion with a spectrosco pic sample of 10^ galaxies as in Model 
I of lMa et aP (l2006l) . The weak lensing constraint on Wp is 
<^(wp) = 0.013, weaker than that from SNela in this idealistic 
scenario, but the area of the 95% contour for weak lensing 
shear tomography is about 30% smaller than that from SNela. 
Aside from serving as a reference, this also demonstrates in 
a tangible way that the criteria of minimizing a(wp) and A 
can result in notably different conclusions. In what remains, 
we make an effort to address constraints in a variety of more 
plausible scenarios. 

The right panel of Fig. |2] depicts Icr constraint contours in 
the more realistic case of imperfect knowledge of the photo- 
metric redshift parameters. As described in §|2] we quantify 
the knowledge of the photometric redshift parameters in terms 
of the size of the spectroscopic sample used to place priors on 
these parameters, A^spec- For the specific results in the right 
panel of Fig.|2l we have chosen A^spec = 3000 and assumed the 
spectroscopic calibration sample to be distributed uniformly 
in redshift over the ra nge 0.1 < z < 1 -7, as might result from 
a forthcoming JDEM (iKim et allllOOl I Albrecht etaD l200g). 
Note that the axes on the right panel of Fig.|2]span a different 
range from the left panel, but that the inner contour showing 
projected constraints from galaxy shear tomography are iden- 
tical in each panel. 

The right panel of Fig. |2] shows a considerable degradation 
of dark energy constraints from either the standard test of lu- 
minosity distance or the dispersion of SNela distance moduli 
alone relative to the case of perfect knowledge of the pho- 
tometric redshift distribution. This demonstrates the impor- 
tance of both a realistic assessment of photometric redshift 
uncertainties and prior information to constrain photometric 
redshifts. However, this panel also depicts the complemen- 
tary nature of these constraints. The luminosity distance test 
alone gives an error on the pivot equation of state parameter of 
(j{wp) w 0.09. Utilizing the information contained within the 
dispersion of SNela as well as the mean decreases the error on 
the pivot equation of state parameter to roughly cr(wp) ~ 0.06 
and decreases the area of the 95% contours relative to ei- 
ther test individually by nearly a factor of ten by breaking the 
prominent degeneracy in the wo-Wa parameter space. Worthy 
of note is the fact that at this level, the constraints from the 
large, photometric SNela sample are comparable to the con- 
straints expected from the spectrosco pic sample used for p ho- 
tometric redshift calibration (see also lAlbrecht et alJ2006l) . so 
they are both interesting and relevant, but the availability of a 
significantly larger spectroscopic sample would obviously be 
more valuable because of the constraints it can place on dark 
energy directly, rather than its ability to calibrate the redshifts 
of a photometric sample. 

Though Fig. |2] shows contours in two simple models, the 
basic point is that SNela dispersions derived from forthcom- 
ing large, photometric surveys are potentially useful. Gener- 
ally, the luminosity distance test is subject to a strong degen- 
eracy between wq, Wa, and JIm- Our choice of priors mitigates 
the influence of J^m, but the wo-Wa degeneracy remains be- 
cause the sample has only a relatively small redshift span and 
therefore limited leverage with which to measure any time 
variation of the dark energy equation of state. Some contem- 
porary studies and several future proposals for SNela-based 
dark energy experiments (such as a Joint Dark Energy Mis- 
sion like SNAP) aim to b reak this de generacy by observing 
high-redshift survey (e.g., Riess et al. 2007). On the contrary, 
our assumed photometric survey has no high-redshift compo- 
nent. In fact, the number of SNela in our assumed survey de- 
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Fig. 3. — Influence of spectroscopic calibration sets on dark energy constraints from a large sample of SNela from a photometric survey. This plot depicts the 
area of the 95% contour in the »'o-»'a plane A, as a function of the size of the spectroscopic SNela sample used to calibrate the photometric redshift model, under 
several assumptions. There are nine lines in the left panel. The solid lines represent constraints as a function of spectroscopic sample size from the luminosity 
distance-redshift test only. The dotted lines represent constraints from the dispersion among SNela distance moduli. The dashed lines show the total constraint 
from both the mean and dispersion in SNela distance moduli. For each line type, the thickest lines show results in a standard case where the spectroscopic 
caHbration sample traces the redshift distribution of SNela in the photometric survey. Alternatively, the thinnest lines show A under the assumption that the 
calibration sample is distributed uniformly in redshift over the interval 7, = 0.l - 1.7. The lines of intermediate thickness show constraints with a spectroscopic 
sample of size Wspcc in addition to a low-redshift "nearby" spectroscopic sample of 500 SNela distributed uniformly in the interval z = 0.03-0.08. The left 
panel shows the case of a known and fixed intrinsic dispersion and the right panel shows the case of a variable intrinsic dispersion. Note that in the case of a 
known intrinsic dispersion, the low-redshift, nearby SNe add little information to the dispersion only test, so the dispersion-only lines in these cases are nearly 
coincident. For clarity, in the right panel, we show only two cases, that of the standard sample and that of the standard sample with an additional, nearby 
spectroscopic sample. 



clines rapidly for 0.6 and only a few percent of the SNela 

in the survey have 0.8 (see lZhan et alj|2008l) . Within the 
photometric survey, the wo-Wa degeneracy is broken by the 
complementarity between the dispersion due to gravitational 
lensing and the standard luminosity distance test. This is a 
potential aspect of complementarity between spectroscopic 
SNela samples and samples from a large, ground-based, pho- 
tometric survey. The latter, bereft of high-redshift (z^ 0.8) 
objects, will be sensitive to w{a) at lower redshifts than pro- 
posed spectroscopic SNela projects that typically extend to 
1.5. Moreover, the dispersion test probes the inhomo- 
geneities in the universe. As such, the dispersion test probes a 
fundamentally different manifestation of dark energy, namely 
its influence on cosmic structure growth, that is not probed by 
the luminosity distance test. 

A further benefit of the additional dispersion information is 
that it helps to break degeneracies between cosmological and 
photometric redshift parameters. In fact, when only consider- 
ing the canonical luminosity distance test, the only constraints 
on the photometric redshift parameters come from the priors 
determined by the photometric redshift calibration set. When 
both pieces of information are exploited, the data provide for 
weak "self-calibration," determining many of the photomet- 
ric redshift parameters with uncertainties a factor of ^2-5 
more stringently than the priors alone. Moreover, Fig. |2] also 
shows that the dispersion test is less sensitive to photometric 
redshift uncertainties and degrades more slowly with decreas- 
ing knowledge of the photometric redshift distribution. This 
useful and interesting result stems primarily from the fact that 
the lensing kernel [W(z',z) in Eq. ©J is inherently broad and 



the dispersion is a comparably mild function of redshift. 

As some spectroscopic calibration of photometric redshifts 
is necessary in order to utilize a photometric SNela sample, it 
is interesting to consider potential constraints both as a func- 
tion of the uncertainty in photometric redshift parameters and 
under a variety of assumptions about what additional SNela 
data will be available. This permits an estimate of the utility 
of SNela from a future large-scale, photometric survey given 
specific photometric redshift calibration programs and addi- 
tional calibration data sets. For this reason, the final two plots 
of this section present results for dark energy parameter con- 
straints as a function of spectroscopic SNela calibration sam- 
ple size A^spec, under a variety of assumptions about additional 
available data. 

First, we consider the area of the 95% ellipse in the wq- 
Wa plane A, as a function of spectroscopic calibration sample 
size, A^spec- This quantity is shown in Figure[3]for several dif- 
ferent combinations of data and model assumptions. In the 
left panel of Fig. [3] we display results for models in which 
the intrinsic dispersion is set to a fixed value and in the right 
panel, for models in which the intrinsic dispersion level is free 
to float. The thickest lines in Fig. [3] show our standard model 
constraints with no nearby SNela sample and a calibration set 
distributed in redshift in a manner identical to the SNela in the 
photometric sample. The intermediate thickness lines show A 
as a function of A^spec in a model that includes the nearby sam- 
ple of SNela. The thinnest lines in Fig. [3]represent the stan- 
dard model with a calibration set that is distributed uniformly 
in redshift in the interval 0.1 < z < 1.7. For illustrative pur- 
poses, we show constraints for spectroscopic sample sizes up 
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Fig. 4. — Constraints on the dark energy pivot equation of state from a large sample of SNela from a pliotometric survey as a function of tlie size of the 
spectroscopically-observed SNela set available to calibrate photometric redshifts. This plot depicts constraint on Wp as a function of the size of the spectroscopic 
sample used Wspec, under several assumptions. There are nine lines in the left panel. The solid lines represent constraints as a function of spectroscopic sample 
size from the luminosity distance-redshift test only. The dotted lines represent constraints from the dispersion among SNela distance moduli. The dashed lines 
show the total constraint from both the mean and dispersion in SNela distance moduli. For each line type, the thickest lines show the results from a standard 
case where the redshift distribution of the spectroscopic calibration sample traces the redshift distribution of SNela in the imaging survey. The thinnest hues 
show (t(»'p) assuming that the calibration sample is distributed uniformly in redshift across the interval 7, = 0.l - 1.7. The lines of intermediate thickness show 
constraints with a spectroscopic sample of size A'spcc in addition to a low-redshift "nearby" spectroscopic sample of 500 SNela distributed uniformly over the 
interval z = 0.03 —0.08. The left panel shows the case of a known and fixed intrinsic dispersion and the right panel shows the case of a variable intrinsic dispersion. 
Note that in the case of a known intrinsic dispersion, the low-redshift, nearby SNe add little information to the dispersion only test, so the dispersion-only lines 
in these cases are nearly coincident. For clarity, in the right panel, we show only two cases, that of the standard sample and that of the standard sample with an 
additional, nearby spectroscopic sample. 



to A^spec =10^; however, we note that such high A^spec are of 
limited practical interest because an upper limit on achievable 
spectroscopic sample sizes is A^spec ^10"* and at such large 
A^spec constraints from the spectroscopic sample itself begin to 
dominate any added information from the photometric sam- 
ple. Note that in Fig. [3] the lines that delineate constraints 
from the dispersion test only in the standard case and in the 
case with an additional nearby SNela sample are nearly coin- 
cident, and may not be readily discernible. The reason for this 
is simple. Under the assumption that the dispersion is known, 
the nearby sample does not serve to calibrate the intrinsic dis- 
persion and adds little information. 

The complementarity of the information contained in the 
mean, /Ii, and the dispersion, cr^.i, is apparent in Fig. [3] The 
information contained in SNela dispersions is not negligible 
and in some cases the constraints from SNela dispersions 
can exceed those from the luminosity distance test alone in 
the regime A^spec^ lO'*. The most interesting cases are likely 
those where the intrinsic dispersion is not known at the out- 
set and/or where there is a nearby SNela sample available. 
In these cases, the information from the standard luminosity 
distance test alone yields a significantly lower area A, than 
the information from the dispersion among distance moduli. 
Combining constraints from both the mean and the dispersion 
information drives A down by roughly a factor of ten com- 
pared to either alone in all cases. This is a manifestation of 
the fact that these two observables are most effective probes 
at different redshifts and so exhibit different degeneracies in 
the wo-Wa plane (see Fig. |2|. In addition. Fig. [3] illustrates 
that constraints from cr^.i degrade more slowly with decreas- 



ing A^spec than do constraints from the standard luminosity dis- 
tance test, a feature already evident in Fig.|2] 

In general, including the low-redshif t, nearby SNe sample 
drives constraints considerably lower dAldering et al.l I2OO2I : 
ICopin et al.l2006tlAlbrecht et alJ200 6'). greatly increasing the 
utility of the luminosity distance test while only yielding a 
slight improvement on the dispersion test i/the intrinsic dis- 
persion level is assumed to be known. Assuming a cali- 
bration set of spectroscopically-observed SNela that is dis- 
tributed uniformly in redshift is moderately more powerful 
than the calibration set that follows the redshift distribution 
of the SNela in the photometric sample. This is because the 
uniform redshift distribution allows for more accurate photo- 
metric redshift calibration at high Zp. This is important for 
two reasons. First, deviations between observables predicted 
by different models are larger at high redshift. Second, the 
redshift distribution of SNela dn/dz is a steep function of 
redshift, dropping as d«/dz cx exp[-32(z-0.5)^] at z 0.5, 
so the scatter of SNela in this regime to higher redshift due 
to redshift errors can lead to significant contamination and 
a Malmquist-type bias. In all cases, A continues to decline 
rapidly with A^spec- In principle, this rate of decline levels off 
as knowledge of photometric redshift distributions become so 
accurate as to allow for the recovery of all available cosmo- 
logical information. This drives all of the constraints to begin 
to converge at large A^spec- In practice, this saturation is not 
achieved until A^spec ~ 3 x 10^, a limit that is so large as to 
be uninteresting (it is larger than the photometric SNela sam- 
pie!). 

The right panel differs from the left in that the intrin- 
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sic dispersion level (Tint, is treated as unknown and permit- 
ted to vary along with A, the cosmological parameters, and 
the photometric redshift parameters. This leads to consid- 
erable degradation in the utility of the dispersion alone to 
constrain dark energy, but in this case the dispersion infor- 
mation constrains the intrinsic dispersion, determining it to 
(T((Tint) ~ 1.6 X 10"^. For comparison, we have evaluated the 
ability of either the nearby sample or a spectroscopic sam- 
ple with with 3000 SNela spread uniformly over the redshift 
interval 0.1 < z < 1.7 and redshift errors of cr-. = lO"'' to con- 
strain (Tint. We find that determinations from such data are 
similar to, but slightly less restrictive, than that from the pho- 
tometric survey after marginalizing over all other parameters. 
Conversely, much of the loss of information due to a varying 
level of intrinsic dispersion is largely recovered in the case 
where the nearby SNela sample is included. Notice that com- 
plementarity between the information in Jli and a^ i leads to 
considerable improvements in A in all cases. 

Another useful measure of the constraining power of any 
set of observables is the error on the equation of state param- 
eter at the pivot redshift, (T(wp). We show the dependence of 
o'(wp) upon A^spec in FigureH) First, note that aiw^) varies over 
a smaller dynamic range than A, so Fig. |4] has an ordinal axis 
with a linear scale as opposed to the logarithmic scale used in 
Fig. El After accounting for this difference there are yet some 
qualitative differences in the behavior of o'(wp) relative to A. 
In the limit of small A^spec and in the absence of a large, low- 
redshift, nearby sample of SNela the dispersion information 
can give useful, independent constraints on Wp. Moreover, 
the dispersion information is generally less useful for improv- 
ing cr(wp) than it is for improving the A figure of merit, par- 
ticularly when the photometric redshift parameters are very 
well known or the nearby SNela sample is available. In these 
cases, A improves because the error ellipses or oriented at a 
significant angle relative to each other, but the shortest axis 
(the constraint on Wp) is largely determined by the luminosity 
distance test alone. Lastly, it is worthwhile noting that in the 
case where there is an available nearby spectroscopic sample, 
the marginal improvement in Wp upon adding additional ob- 
jects to the spectroscopic photometric calibration sample is 
relatively small beyond A^spec ^ 2 x 10"^. 

To this point, we have considered cases in which the intrin- 
sic dispersion of SNela does not vary with redshift. It is not 
unlikely that the intrinsic dispersions in SNela will change 
with redshift. For example, SNela will be observed in dif- 
ferent bands relative to their rest frames at different redshifts 
with the consequence that st andard candle calibration will be 
a fun ction of redshift (e.g., iHamuy et al.llI995l: iPrieto et aP 
|2006|) . In addition, there is evidence that the la class of su- 
pernovae is composed of distinct sub-classes and the rela- 
tive m ix of these sub-classes is expected to change with red - 
shift ( Hamuv etal.lll995t lHowellll200lt iHowell et al.ll2007h . 
Clearly, if the intrinsic dispersion is allowed to be an arbi- 
trary function of redshift, the utility of the added dispersion 
information is completely eliminated as a cosmological probe. 
However, we have considered a less pathological case of a 
redshift-dependent intrinsic dispersion that is monotonic and 
follows a power-law. 

To be specific, we considered a model with CTint(z) = 
o'int.o(l with both (Jint.o and (3 parameters that are fixed 
by the data (/3 = in the fiducial model about which we per- 
turb). We do not show the results in Fig.|3]or Fig.|4]for clarity, 
but note that further degradation is not devastating. In particu- 
lar, in the case with no nearby sample, further degradation by 



allowing for crint(z) to vary as a power law in (1 + z) is a factor 
of ~ 2.5 in and ~ 18% in (7(wp) relative to the case of con- 
stant (Tint. With a nearby spectroscopic sample, (Jiatiz) at low 
redshift can be effectively determined and the degradation rel- 
ative to the case of a redshift-independent intrinsic dispersion 
is only ~ 60% in A and just ^ 10% in cr(vvp). In addition, 
these data would limit (3 with uncertainties of cr(/3) ~ 0.22 
and a((5) ~ 0.09 in cases with and without the nearby sam- 
ple respectively, a level that is again comparable to the fore- 
casts from a spectroscopic sample. Moreover, it is clear 
that the lensing contribution to the dispersion must be ac- 
counted for in order to extract the intrinsic SNela properties. 
In summary, though the information loss in cases of redshift- 
dependent intrinsic dispersion is not insignificant, SNela in 
a large, photometric survey can bring significant constrain- 
ing power to bear on dark energy even with moderate, and 
poorly-understood time variation of the intrinsic dispersion. 
Of course, to the degree that the distribution of distance mod- 
uli varies markedly and rapidly with redshift outside of the 
expected range, this distribution will provide a useful han- 
dle on the SNela properties themselves, including perhaps 
the e volution of SNela and SNela population demographics 
(e.g.JHamuv et al.l[T99llHowellll200UlM annucci et al."2006t 
ISuUivan et al.ll2006t iHoweUetalTEoOTL lSarkar et al. 2008a). 

4. IMPLICATIONS AND CAVEATS 

We have studied simultaneous constraints on dark energy 
coming from both the evolution of the mean luminosity dis- 
tance of SNela as a function of redshift and the dispersion 
among SNela fluxes in a SNela data set that may arise from a 
large, wide, and fast photometric survey such as that proposed 
for the LSST. Sources of dispersion among the measured 
fluxes of SNela in such a survey include intrinsic disperision 
among SNela (including standard candle calibration), uncer- 
tainty in photometric redshifts, and magnification due to grav- 
itational lensing. We have shown that the additional disper- 
sion information complements the traditional luminosity dis- 
tance test in several ways. 

First, the dispersion information breaks a degeneracy be- 
tween the contemporary dark energy equation of state wq, and 
dark energy equation of state evolution, parameterized here 
by the common, benchmark parameter Wa. Using the lumi- 
nosity distance (mean flux) test alone, this degeneracy can be 
broken by observing high-redshift SNela to increase the lever 
arm of the data in redshift. The dispersion information leads 
to constraints that break the degeneracy between wq and Wa 
without the nead for a high-redshift (z^ 0.8) SNela sample 
(see Fig.|2]l. 

Second, photometric surveys require some calibration of 
photometric redshifts. Including dispersion information may 
allow for mild internal self calibration of uncertainties in the 
true redshift distribution of observed SNela given their pho- 
tometric redshifts. In the context of our models, this results 
in constraints on photometric redshift model parameters (see 
§12]) that are a factor of ~ 2-5 more constraining than the pri- 
ors on these parameters from a spectroscopic sample of size 
A^spec = 3000. The realized level of improved calibration de- 
pends upon the details of both the model and the additional 
SNela calibration data that are available, but in all cases the 
luminosity distance test alone does not serve to calibrate the 
photometric redshift distribution. 

Lastly, it is entirely possible that there will be evolution in 
the distribution of SNela intrinsic luminosities that signifi- 
cantly degrade the forecasts for dark energy constraints that 
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we present (see § [3]). All approaches to SNela cosmology, 
even those that exploit the traditional luminosity distance test 
alone, rely upon some knowledge and/or assumptions about 
the intrinsic distribution of SNela luminosities to obtain cos- 
mological constraints. In cases where the evolution of the 
spread in luminosities at fixed redshift is significant over the 
redshift range 0.2< 0.8, a photometric survey cannot use 
lensing to constrain dark energy directly. However, we have 
shown that such a survey can measure the evolution of the in- 
trinsic SNela luminosity distribution at levels comparable to, 
or better than, a spectroscopic survey containing a few thou- 
sand SNela. This should aid the luminosity distance test by 
providing a better understanding of possible SNela evolution 
and demographics and, therefore, a better understanding of 
potential parameter biases and uncertainties. 

As we have mentioned in the preceding paragraph, in exam- 
ining a photometric SNela sample, it is necessary to consider 
uncertainties in the photometric redshift distribution. We have 
assumed that some calibration will be done with a spectro- 
scopic sample of an unspecified size, A^spec- Figure[3]gives the 
expected area of the 95% ellipse in the wq-w-^ plane. A, as a 
function of A^spec- A decreases rapidly with A^spec in all cases 
because increasing A^spec allows access to information over a 
broader range of observed redshifts and helps to break degen- 
eracies in the wq-w^ plane. Figure |4] shows the constraints on 
the equation of state parameter at that redshift where the data 
best constrain it, the so-called pivot equation of state Wp, as 
a function of A^spec- This figure shows that the marginal im- 
provement in (T(wp) with increasing A^spec decreases beyond 
■^spec ^ a few X lO-'. This result is interesting because the up- 
per limit on the achievable size of a spectroscopic sample is 
near A^spec ~ lO"* and, moreover, with A^spec^ lO'* constraints 
from the spectroscopic sample alone (which has miniscule 
redshift errors, ^ 10""*) would begin to dominate over that 
achievable with a large, photometric sample even in the ideal 
limit where statistics, and not systematics, dominate the error 
budget. Both Fig. |3]and Fig. |4] illustrate that at fixed A^spec 
it is more fruiful to employ a spectroscopic calibration set 
weighted more toward high-redshift (z^ 0.5) SNela than the 
distribution of SNela from the imaging survey. 

The previous paragraph closed with an allusion to system- 
atic error, and indeed sources of systematic error are a major 
caveat to the present study. The main advantage of a photo- 
metric sample is its size and the additional statistical leverage 
that come with a large sample size. However, once system- 
atic errors dominate the error budget, it is no longer possi- 
ble to take advantage of this benefit. What is more, a pho- 
tometric survey places emphasis on many potential sources 
of systematic error. We have accounted for some of these 
complications. For example, we have allowed for consistent 
determination of the photometric redshift distribution and for 
variation of the average SNela luminosity and intrinsic dis- 
persion, including mild evolution. This is indicative that our 
results will be somewhat robust to evolution in the intrin- 
sic SNela population (e.s. iHamuy et aTI|1995t lHowellll2001t 
ISuUivan et al.1l2006t iHowell et al.ll2007t TSarkar et al.ll2008al) . 



Other sources of systematic error are likely to be important. In 
fact, the SNAP collaboration, which aims to observe and ex- 
ploit a spectroscopically-observed SNela sample from space, 
has assumed a fundamental limit in the ability to calibrate 
supernovae to within Sfi ^ 0.02 mag over a reds hift differ- 
ence of & ~ 0.1 to specify their survey (e.g. lRim~e t al. 200'^. 
Yet it remains to be seen whether such a systematic limit will 



be realized, and it is possible that it may be higher or lower. 
Further, without spectroscopy, not only are redshifts uncer- 
tain, but supernovae type identification, extinction and K- 
corrections, a nd standard candle cahbration all be come more 
difficult (e.g. 'Pinto et al."2004| [MetoeLaLl 11006'). We have 
not explicitly accounted for any of these, largely because sys- 
tematic levels are uncertain and modeling them is both com- 
plex and dependent upon survey strategy. These complica- 
tions place a detailed treatment of systematics beyond the 
scope of our study (the studies by the DETF, as well as those 
of Zhan et al. 2008 and lHannestad et al . 2008 neglected a de- 
tailed treatment of systematics as well). If the assumed SNAP 
systematics floor is realized, our forecasts for dark energy 
equation of state constraints increase by a factor of ^ 3, sig- 
nificantly reducing the additional cosmological information 
that may be extracted from the photometric sample of SNela. 

As described in §|2] the calculations we have performed to 
estimate the additional information contained in the variety of 
SNela distance moduli at a fixed redshift are greatly simpli- 
fied. Our idealization may lead to an over-estimate of the con- 
straining power of SNela. However, we have also neglected 
some potentially-available information. In particular, the dis- 
tribution of lensing magnifications is highl y non-Gaussian 
dSchnei der & Wagoner 1987; Sasaki 1987; W ambsganss et al.l 
119971; iMunshi & Jainii2000c iWang et al.. 2002) and we have 
neglected any information beyond the dispersion because in- 
cluding the full information is computationally-intensive and 
current theoretical estimates of the magnification distribution 
are inadequate for our application. The left panel of Fig- 
ure|5]shows the non-Gaussianity of the lensing distribution at 
a fixed redshift z = 0.8 using the fitting formula of Wang et ai] 
( |2002|) . The right panel of Figure |5] shows the distribution of 
distance moduli in a bin of width 6z = 0.03 centered at red- 
shift z = 0.8. In our fiducial model, the bin contains roughly 
n w 19,000 SNela and Fig. |5] demonstrates that the non- 
Gaussianity of the lensing distribution would be detectable 
even by examining the distribution of distance moduli in this 
single bin. For the realization shown, the non-Gaussianity 
of the lensing distribution is detectable at slightly more than 
^ 2(7 from this bin alone. This strongly suggests that the 
lensing deviation from Gaussianity may be detectable in a 
complete analysis of a future photometric survey, where the 
full shape of the magnification distribution would be utilized 
and where all supernovae could be considered (which requires 
better predictions). This is important for two reasons. First, 
the unique distribution of fluxes induced by lensing may aid 
in distinguishing the gravitational lensing effect from other 
additional sources of dispersion and from evolution in the in- 
trinsic properties of SNela. Second, to the degree that this 
non-Gaussianity can be detected, it will bring additional in- 
formation to bear and improve dark energy constraints. 

To conclude, a future imaging survey that scans a large frac- 
tion of the sky rapidly will enable the discovery of ~ 10^-10^ 
SNela. The survey of the LSST is the canonical exam- 
ple of such an endeavor The number of SNela with light 
curves of sufficient quality to be used for cosmology depend 
upon specific survey strategies; however, such a SNela sur- 
vey may be able to constrain dark energy properties using 
both the traditional luminosity distance test and the spread 
of supernovae apparent brightnesses as a function of photo- 
metric redshift. Unfortunately, current theoretical treatments 
of gravitational lensing, in particular estimates of the magni- 
fication distribution, are not reliable enough to analyze any 
such data set. In the coming years, it will be necessary to 
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Fig. 5. — The non-Gaussianity of the leasing distribution. Both panels show the distribution of the shift in distance modulus of SNela. The left panel 
is illustrative. The solid, smooth curve shows the distribution of distance moduli for sources at redshift z = 0.8 due only to lensing. This result was computed 
using the fitting formula of Wang et al. (2002). We display this quantity at relatively high redshifts only because the formula of Wang et al. (2002) is no longer 
internally self-consistent below 7, ^ 0.6. The dotted curve shows a Gaussian distribution with the same mean and dispersion as the lensing distribution. The 
Gaussian distribution is truncated at the minimum magnification, con'esponding to a maximum shift in distance modulus 5^niax. The histogram shows the 
distribution of distance modulus shifts from a random realization of 19,000 SNela at z = 0.8. The right panel shows something closer to what would be seen 
in any observed set of SNela. In this panel, we show the distribution of distance moduli (relative to the mean) for a random sample of SNela in a bin of width 
5z = 0.03 at z = 0.8. The sample consists of 19, 000 SNela, which is roughly the number of SNela expected in such a bin according to our fiducial model described 
in §|2] In this case, the distribution includes the intrinsic dispersion, the dispersion caused by poorly-controlled photometric redshifts, and the lensing distribution. 
For clarity, the probability distribution is plotted relative to the same distribution under the assumption of the lensing contribution is a Gaussian, Poaussian Lcn.sl'^M)- 
The error bars reflect the statistical error on P(5fi) from a sample of this size and portray of the detectability of the unique shape induced by lensing in the absence 
of significant errors. 



refine these estimates, including p otential theoretical uncer- 
tainties dHuterer & Takada 2005; Zhan & Knoxll2004t IWhlt^ 
120041: iRudd et al.i i2008i) . In the end, the utility of such a 
sample to constrain dark energy is subject to limitations of 
systematic error and SNela evolution. However, if observa- 
tional systematics or SNela evolution prove to be limiting, 
such a large photometric survey will provide useful informa- 
tion about these issues. In either case, an LSST-like imaging 
survey will revolutionize SNela cosmology. 
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